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ABSTRACT 

Aims. Having established the binary status of nineteen O-type stars located in four Cygnus OB associations, we now determine 
their fundamental parameters to constrain their properties and their evolutionary status. We also investigate their surface nitrogen 
abundances, which we compare with other results from the literature obtained for galactic O-type stars. 

Methods. Using optical spectra collected for each object in our sample and some UV data from the archives, we apply the CMFGEN 
atmosphere code to determine their main properties. For the binary systems, we have disentangled the components to obtain their 
individual spectra and investigate them as if they were single stars. 

Results. We find that the distances of several presumably single O-type stars seem poorly constrained because their luminosities are 
not in agreement with the ’’standard” luminosities of stars with similar spectral types. The ages of these O-type stars are all less than 
7 Myrs. Therefore, the ages of these stars agree with those, quoted in the literature, of the four associations, except for CygOB8 for 
which the stars seem older than the association itself. However, we point out that the distance of certain stars is debatable relative 
to values found in the literature. The N content of these stars put in perspective with N contents of several other galactic O-type 
stars seems to draw the same five groups as found in the ’’Hunter” diagram for the O and B-type stars in the LMC even though their 
locations are obviously different. We determine mass-loss rates for several objects from the Ho: line and UV spectra. Finally, we 
confirm the ’’mass discrepancy” especially for O stars with masses smaller than 30 M e . 

Key words. Stars: early-type - Stars: binaries: spectroscopic - Open clusters and associations: individual: Cygnus OB 1 - Open clusters 
and associations: individual: Cygnus OB3 - Open clusters and associations: individual: Cygnus OB8 - Open clusters and associations: 
individual: Cygnus OB9 


1. Introduction 

Massive stars are the engines of the galaxies. This small part 
of the stellar population has a considerable influence on its envi¬ 
ronment both mechanically and by radiation. Indeed, their strong 
ionizing fluxes power the H ii regions. Furthermore, their power¬ 
ful winds and their deaths as supernovae enrich their surround¬ 
ings with heavy chemical elements and reshape the interstellar 
medium, triggering new generations of stars. Because their life¬ 
time is relatively short, massive stars are excellent tracers of star 
formation. Therefore, constraining their properties gives us the 
opportunity to better understand how they form and how they 
evolve. 

The fates of these stars are mainly governed by their initial 
mass, but also to some extent by their rotation and their mass 
loss. In addition to giving more oblate shapes to the stars, the 
rotation modifies their temperature gradient at the surface and 
also induces an internal mixing that affects the angular momen¬ 
tum and chemical element transport. The rotation thus influences 
the surface abundances and the lifetimes of the massive stars 
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on the main sequence dMevnet & Maedeif200(j) . Several obser¬ 
vational st udies analysed thes e two aspects of rotation. On the 
one hand, [Hunter et al.1 62 0081120091) reported that a large part 
of B-type stars in the Large Magellanic Cloud (LMC) presented 
a clear trend between their projected rotational velocity (v sin i) 
and their surface nitrogen content. This trend was detected in 
about 60% of the B star population. The remaining 40% consists 
of B stars with slow rotation and high nitrogen enrichment and 
unenriched B stars with fast rotation. These groups are currently 
not explained by the evolutionary models of single massive stars 
which include rotational mixing. On the other hand, numer¬ 
ous investigations have shown that the rotation brings modifi¬ 
cations to the evolu tionary tracks and to the isochrones of mas¬ 
sive st ars (see e.g. iMevnet & Maedetl l2000l 12003b iBrott et ail 
1201 lal) . Besides rotation, the stellar winds also play a crucial 
role in the evolution of massive stars. Indeed, these objects lose 
a significant fraction of their mass during their life. However, 
the exact quantity of ejected mass by a massive star remains 
difficult to determine because the true structure of these winds, 
and no tably their homogeneit y, is still poorl y known. Although 
direct (lEversberg et al.l fl998) and indirect feouret et al.l [2005 ) 
evidence highlighted the presence of outward-moving inhomo¬ 
geneities in the winds of O-type stars, the shapes, the sizes, and 
the optical depths of these inhomogeneities (or clumps) still re- 
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main unanswered questions. Investigations of the structure of 
these winds seem, however, to indicate that the clumps are spher¬ 
ical (ISundavist et alJl201 It Herve et al.ll 2012l) rather than flat¬ 
tened as proposed bv iFeldmeier et al! ~d2003l) . The presence of 
clumps in the stellar winds thus modifies the determination of 
mass-loss rates of massive stars since many mass-loss diagnos¬ 
tics are dependent on the density squared (such as notably Ha). 


To bring more constraints on these properties, we have de¬ 
termined and analysed the fundamental parameters of nineteen 
O-type stars belonging to four OB associations in the Cygnus 
complex. This area is an active star-forming region that includes 
a huge number of massive stars. These targets have already been 
presented in iMahv et ah ( 2013L hereafter Paper I) who estab¬ 
lished their multiplicity and analysed the distribution of the or¬ 
bital parameters of the binary systems. Because we know their 
multiplicity, in the present paper we can be more accurate on the 
determination of valuable physical parameters of these objects. 
Although all these targets are formed in the Cygnus region, they 
constitute a non-homogeneous sample of main-sequence band, 
giant, and supergiant massive stars. Among this sample of stars, 
we also detected in Paper I four binary systems: three SB2s and 
one SB1, all with orbital periods shorter than 10 days. Therefore, 
for the SB2s, we disentangled the observed spectra of these sys¬ 
tems to obtain the individual spectrum of each component to in¬ 
vestigate it with an atmosphere code as we do for single stars. 


In the following section, we present the observations of these 
targets. The atmosphere code and the main UV/optical diagnos¬ 
tics are described in Sect. 0 The results and a discussion are 
given in Sect.[4]and in Sect. [3 respectively. Finally, our conclu¬ 
sions are provided in Sect. [6] 


2. Observations and data reduction 

We focus on the same data as presented in Paper I. In addition, 
we retrieve the IUE spectra (SWP with a dispersion of 0.2 A) 
to model the stellar winds of five objects. These data were ob¬ 
tained between 1979 and 1992 with exposure times from 650 s 
to 8400 s. These spectra cover the [1200-1800] A wavelength 
domain, containing the C iv 1548-1550 and N v 1240 resonance 
P-Cygni profiles as well as the Ov 1371 line, used to estimate 
the wind clumping in early O stars. We also add to our sam¬ 
ple four FUSE spectra retrieved from the archives. Because of 
the interstellar absorption, we only focus on the [1100-1200] A 
wavelength region, and mainly on the P v 1118-1128, C iv 1169, 
and C hi 1176 lines. Unfortunately, the FUV and UV spectra are 
not available for all the stars in our sample. Therefore, the mass- 
loss rates of the stars for which none of these spectra exists will 
only be derived on the basis of the Ha line. 

For the targets detected as binary systems in Paper I, we 
used the disentangled spectra corrected for the brightness ratio 
which were already prese nted in Paper I. We note th at we used 
a programme based on the Gonzalez & Levato ( 200 61 technique 
to compute the individual spectra of both components of a bi¬ 
nary system, but also to refine the radial velocities by applying 
a cross-correlation technique. These spectra can thus be consid¬ 
ered as mean spectra for the components of the binary systems. 
To have as many diagnostic lines as possible, we favour the op¬ 
tical spectra with the widest wavelength coverage, as well as the 
highest spectral resolution and the highest signal-to-noise ratio. 
The journal of optical and FUV/UV observations is presented in 
Table □ 


3. Modelling 

We use the code CMFGEN dHillier & Miller! 1998t) for the quan¬ 
titative analysis of the optical and FUV/UV spectra. CMFGEN 
provides non-LTE atmosphere models including winds and line¬ 
blanketing. CMFGEN needs as input an estimate of the hydro- 
dynamical structure that we construct from TLUSTY models 
(taken from the OSTAR2002 grid of lLanz & Hubeny|2003[) con¬ 
nected to a /? velocity law of the form v' = v 00 (l - R/r'f , where 
Voo is the wind terminal velocity. For the stars having only opti¬ 
cal spectra, we adopted/?_= 0.8, which represents a typical value 
for O dwarfs (see e.g. lRepolust et al.l2004l) . Our final models in¬ 
clude the following chemical elements: Hi, Hei-n, Cii-iv, Nn- 
v. On-vi, Neii-ni, Mgn, Sin-iv, Siii-vi, Piv-v, Arm-iv, Aim, 
Fe ii- vii, and Ni iii-v with the solar composition of iGrevesse et al.l 
( 2007 1 unless otherwise stated. CMFGEN also uses the super¬ 
level approach to reduce the memory requirements. On average, 
we include about 1600 super levels for a total of 8000 levels. 
For the formal solution of the radiative transfer equation lead¬ 
ing to the emergent spectrum, a microturbulent velocity varying 
linearly with velocity from 10 km s 1 to 0.1 x Voo was used. We 
include X-ray emission in the wind since this can affect the ion¬ 
ization balance and the strength of key UV diagnostic lines. In 
practice, we adopt a temperature of three million degrees and 
we adjust the flux level so that the X-ray flux coming out of 
the atmosphere matches the observed Ux/Tboi ratio. We sim¬ 
ply adopt the canonical value Lx/Tboi = 10 7 (ISana et al.ll2006b 
lNazel2009h . In practice, we proceed as follows to derive the stel¬ 
lar and wind parameters. 

- Effective temperature: we use the classical ratio between the 
strengths of He i and He ii lines to determine T e ff. The main 
indicators are the He i 4471 and He n 4542 lines, but addi¬ 
tional diagnostics can be built with the He i 4026, He i 4389, 
Hei 4713, Hei 4921, Hei 5876, Hen 4200, and Hen 5412 
lines. When possible we also use the C iv 1169 to C m 1176 
line ratio, which has_been shown to provide a useful temper¬ 
ature diagnostic dHeap et al.l 12006 1. The typical uncertainty 
on the r eff determination is 1000 K, except for the binary 
components for which the uncertainty on 7/n' is 2500 K. 

- Gravity: the wings of the Balmer lines H/j, Hy, and He) are 
the main indicators of log g. Generally, an accuracy of about 
0.1 dex on log g is achieved. However, because of side effects 
of the disentangling programme, log g of the binary compo¬ 
nents has an uncertainty of 0.25 dex. 

- Wind terminal velocity: the blueward extension of UV P- 
Cygni profiles provides Voo+v max where v max is the maximum 
turbulent velocity. Using the above relation for microturbu¬ 
lent velocity gives us a direct determination of with an 
accuracy of 100 km s 1 . For stars without UV spectra, we 
use as input the standard values given bv IPrinia et alJ ( 1990 1 
according to their spectral types. 

- Mass-loss rate: we use two diagnostics to constrain the mass- 
loss rate, the UV P-Cygni profiles and Ha. For the UV do¬ 
main, we focus mainly on Nv 1240, Siiv 1394-1403, and 
Civ 1548-1551, and when FUSE spectra are available, we 
consider C iv 1169 and C hi 1176 as additional diagnostics to 
those lines for mass-loss rate determination. A single value 
of M should allow a_good fit of both types of lines even 
though a study by iMartins et all ( 2012 ) has shown that this 
was not always the case. In our analysis, we reach a general 
agreement between UV and Ha lines by playing on other 
parameters such as the clumping or the /?. However, when 
no FUV/UV data are available, we fix both values to 1.0 
(homogeneous model) and 0.8 fsee lRenolust et all2004l) . re- 


2 











































Table 1 . Journal of optical and UV observations 



Optical 




IUE 



FUSE 


Star 

Instrument Spectral range 
[A] 

HJD 

id] 

Obs. date 

data ID 

HJD 

[d] 

Obs. date 

data ID 

HJD 

Obs. date 


Cyg OBI 


HD 193443 

- 

[4450-4900] 

disentangled spectra 

- 

- 

- 

- 

- 

- 

HD 228989 

- 

[4450-4900] 

disentangled spectra 

- 

- 

- 

- 

- 

- 

HD 229234 

Espresso 

[4000-6750] 

2454 992.7849 

10 Jun 2009 

- 

- 

- 

- 

- 

- 

HD 193514 

Elodie 

[4000-6800] 

2453 600.4087 

17 Aug 2005 

SWP18145 

2445 241.2984 

28 Sep 1982 

E0820701 

2453 165.3275 

08 Jun 2004 

HD 193595 

Espresso 

[4000-6750] 

2454 989.8582 

07 Jun 2009 

- 

- 

- 

- 

- 

- 

HD 193682 

Espresso 

[4000-6750] 

2454 989.9081 

07 Jun 2009 

SWP09022 

2444 376.4267 

16 May 1980 

E0820301 

2453 365.5448 

09 Jun 2004 

HD 194094 

Espresso 

[4000-6750] 

2455 820.7063 

16 Sep 2011 

- 

- 

- 

- 

- 

- 

HD 194280 

Aurelie 

[4450-4900] 

2454 407.2709 

02 Nov 2007 

- 

- 

- 

- 

- 

- 

HD 228841 

Espresso 

[4000-6750] 

2454 989.8092 

07 Jun 2009 

- 

- 

- 

- 

- 

- 

Cyg OB3 

HD 190864 

Elodie 

[4000-6800] 

2452 134.4733 

12 Aug 2001 

SWP06946 

2444168.7353 

22 Oct 1979 

E0820501 

2453 148.7486 

23 May 2004 

HD 227018 

Aurelie 

[4450-4900] 

2454 717.3188 

07 Sep 2008 

- 

- 

- 

- 

- 

- 

HD 227245 

Espresso 

[4000-6750] 

2455 727.7497 

15 Jun 2011 

- 

- 

- 

- 

- 

- 

HD 227757 

Espresso 

[4000-6750] 

2455 726.7373 

14 Jun 2011 

- 

- 

- 

- 

- 

- 

Cyg OB8 

HD 191423 

Elodie 

[4000-6800] 

2453 247.3339 

29 Aug 2004 

SWP16212 

2445 000.2270 

30 Jan 1982 

E0821301 

2453 165.7979 

09 Jun 2004 

HD 191978 

Elodie 

[4000-6800] 

2453 247.4196 

29 Aug 2004 

SWP46532 

2448 976.0517 

19 Dec 1992 

- 

- 

- 

HD 193117 

Espresso 

[4000-6750] 

2454 991.9784 

09 Jun 2009 

- 

- 

- 

- 

- 

- 

Cyg OB9 

HD 194334 

Espresso 

[4000-6750] 

2455 053.7053 

10 Aug 2009 

- 

- 

- 

- 

- 

- 

HD 194649 

- 

[4450-4900] 

disentangled spectra 

- 

- 

- 

- 

- 

- 

HD 195213 

Aurelie 

[4450-4900] 

2454 717.5787 

07 Sep 2008 

- 

- 

- 

- 

- 

- 
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spectively. We also stress that some objects were only ob¬ 
served on [4450-4900] A wavelength domain. Since no suf- 
bcient diagnostic exists to determine the mass-loss rate on 
this region, we also use as input the values provided by 
iMuiires et al] (12012h as mass-loss rate values for our mod¬ 
els. 


Given the complexity of the parameter space for CMFGEN 
models, it is not possible for us to describe all the effects of every 
parameter change o n the synthetic spec t ra, but we refer to dedi¬ 
cated papers su ch as lHillier et al.l (120031) . iBouret et al.l (l2005h . or 
iMartinsI (1201 lh . In our analysis, we decide to vary the different 
parameters until we obtain the solution which provided the best 
X 2 fit. To this end we generated a non-uniform grid composed of 
between 10 and 40 models depending on the wavelength cover¬ 
age of the observed spectra (and thus of the number of param¬ 
eters to determine). From these best-fit models, we reproduce 
the UV-Optical-Infrared spectral energy distribution (SED) for 
each star (Fig.0 from the UBVJHK fluxes (see Table [2} and we 
include the UV fluxes for the stars that have FUSE and/or IUE 
spectra to constrain the distance and the extinction. The galactic 
reddening law of tCardelli et al.l (1989|) is used. We also derive Ry 
when the UV spectr um is available. When this is not the case, we 
use Ry values from Patriarchi et al.l (1200 lL l2003h and if no value 
exists, we fix it to 3.1. We then recompute the My with the new 
extinction and the new distance values and we derive the lumi¬ 
nosities on the basis of the My and of the bolometric corrections 
given in [Martins & Plez ( 2006 ). Once the new luminosity is ob¬ 
tained, we compute a new model with this value to improve the 
best-fit models. The uncertainties calculated on the luminosities 
are mainly due to the poorly-known distances. To quantify these 
uncertainties, we computed the largest difference between the 
obtained luminosity and the luminosity computed with the min¬ 
imum or maximum mean distances of the associations listed in 
[Humphreys! (fl978h . 

The projected rotational velocities of the stars are obtained 
from the Fourier transform method (iSimon-Dfaz & Herrerol 
12007 ). The synthetic spectrum is then convolved by v sin /. 
The rotationally broadened synthetic line profiles usually pro¬ 
vided a good match to the observed profile for the fast rota¬ 
tors. In the case of moderately to slowly rotating objects, some 
amount of macroturbulence has to be introduced to correctly 
reproduce the line profiles of several features (e.g. Hei 4713, 
C iv 5812 , He 1 5876). The need for extra broaden i ng is well doc¬ 
umented dHowart h et al] 19971: iRvans et aljl2002tlHowarth et al 


2007: iNieva & Przv billa 2007; Martin s et al. 2010|; iFraser etal 
2010 ). but its origin is unclear. lAerts et akT i 2009h notably sug¬ 


gested that non-radial pulsations could trigger large-scale mo¬ 
tions (hence macroturbulence), but this needs to be confirmed 
by a study cove ri ng a wider parameter space. A study by 
ISimon-Dfaz et alj d2010t) also showed that the amplitude of 
macroturbulence was correlated to the amplitude of line profile 
variability in a sample of OB supergiants. In practice, we in¬ 
troduce macroturbulence by convolving our rotationally broad¬ 
ened synthetic spectra with a Gaussian profile, thus mimicking 
isotropic turbulence. This is obviously a very simple approach, 
but it significantly improves the quality of the fits. In practice, we 
used He i 4713 as the main indicator of macroturbulence since it 
is present with sufficient signal-to-noise ratio in all the stars in 
our sample. Moreover, secondary indicators are the Oni 5592 
and He i 5876 lines. 

To determine the nitrogen content of our stars, we rely 
mainly on the N in lines between 4500 and 4520 A. They are 
present in absorption in the spectra of all stars, they are not af¬ 


fected by winds and they are strong enough for the abundance 
determination. The uncertainties are of the order of 50%. They 
are estimated from the comparison of the selected N m lines to 
models with various N content. The uncertainties do not take any 
systematics related to atomic data into account. 

When UV spectra are available (both FUSE and IUE), the 
degree of inhomogeneities of the stellar winds is determined. 
Clumping is implemented in CMFGEN by means of a volume 
filling factor / following the law / = /«, + ( 1 — f 00 )e~ v ^ d , where 
foe is the maximum clumping factor at the top of the atmosphere 
and I'd a parameter indicating the position where the wind starts 
to be significantly clumped. As shown by IBouret et al. ( 2005 ). 
Pv 1118-1128, Ov 1371, and Niv 1720 are UV features es¬ 
pecially sensitive to wind inhomogeneities. We use these lines 
to constrain f m in the stars in our sample. Unfortunately, the 
last two lines are mainly good indicators for early O-type stars. 
Therefore, we are not able to use these two lines to constrain the 
clumping factor for HD 191978 given the absence of a FUSE 
spectrum. 


4. Results 

The derived stellar and wind parameters are listed in Tabled The 
CMFGEN best-fit models are given in the Appendix. Below, we 
briefly comment on each star. 

4.1. The CygOBI association 

HD 193443: In Paper I, we estimated the brightness ratio to 
3.9 + 0.4. The disentangled spectra corrected for this ratio are 
relatively well fitted. However, the wavelength coverage does not 
allow us to determine the wind parameters. The individual stellar 
parameters of the primary are more reminiscent of a giant star, 
whilst the secondary has the characteristics of a main-sequence 
object. This confirms the luminosity classes of both components 
mentioned in Paper I. 

HD 228989: We found in Paper I a brightness ratio of 1.2 + 
0.1 between the primary and the secondary. The spectral lines 
of the disentangled spectra are relatively well fitted even though 
the He n 4686 line of both stars stronger than with the synthetic 
profiles. This is probably due to a poor estimation of the wind 
parameters, but, as stressed for HD 193443, the absence of other 
diagnostic lines prevents us from providing a good determination 
of the wind parameters. 

HD 229234: This object was reported as SB1 in Paper I. 
However, since the signature of the putative companion is not 
detectable in the observed spectra, we consider that the these 
spectra are only produced by the primary star, suggesting that 
the derived parameters correspond to an upper limit. The fit is 
relatively good and indicates stellar parameters similar to those 
of a giant star. This agrees with the spectral classification estab¬ 
lished in Paper I. 

HD 193514: The param eters derived for HD 193514 are sim¬ 
ilar to those of lRenolust et al.ld2004l) . However, the wind param¬ 
eters are slightly different. Indeed, we need a clumping factor of 
0.01 to correctly reproduce the FUV lines. This thus affects the 
mass-loss rate of the star even though Mna/ \/7 remains quite 
close in both analyses. The emissions of the N hi 4634^-1 lines 
in the spectrum of HD 193514 are poorly reproduced, but it is 
probably not due to the N content, given that the triplet at 4500- 
4520 A and the N iv 1720 line are well fitted. Moreover, we stress 
that an increase of the mass-loss rate does not improve the qual¬ 
ity of the fit. 
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Table 2. UBVJHK photometric parameters of the stars in our sample 


Star 

Cluster 

U 

B 

V 

J 

H 

K 

R v 

E(B-V) 

Dist. 

(SED) 

My 

Dist. 

Humphreys (1978) 

Cyg OBI 

HD 193443 


7.08° 

7.61 

7.26 

6.40 

6.34 

6.34 

3.10 

0.61 

- 

-5.93 

1.82 

HD 193514 


7.32" 

7.84 

7.43 

6.35 

6.24 

6.18 

3.00 

0.74 

1.71 

-5.96 

1.90 

HD 193595 

Ber86 

8.50* 

9.08 

8.78 

7.92 

7.87 

7.85 

3.02 

0.69 

1.85 

-4.64 

1.82 

HD 193682 


8.39* 

8.87 

8.42 

7.30 

7.19 

7.14 

2.82 

0.83 

1.95 

-5.37 

1.82 

HD 194094 


9.18* 

9.57 

9.05 

7.74 

7.59 

7.58 

2.99 

0.89 

2.65 

-5.73 

2.42 

HD 194280 


8.94° 

9.08 

8.42 

6.74 

6.60 

6.49 

2.95 

1.06 

2.67 

-6.84 

2.20 

HD 228841 

Ber86 

9.08" 

9.50 

9.01 

7.91 

7.78 

7.75 

2.78 

0.88 

2.70 

-5.59 

2.09 

HD 228989 

Ber86 

10.29* 

10.49 

9.83 

8.15 

7.97 

7.90 

3.10 

0.92 

- 

-3.81 

1.43 

HD 229234 

NGC6913 

9.35 d 

9.58 

8.94 

7.34 

7.18 

7.10 

2.87 

1.07 

1.50 

-5.02 

1.59 

CygOB3 

HD 190864 

NGC 6871 

7.20" 

7.93 

7.78 

7.26 

7.28 

7.26 

2.45 

0.58 

2.15 

-5.30 

2.29 

HD 227018 

NGC6871 

8.73* 

9.35 

9.00 

8.04 

7.95 

7.90 

3.28 

0.71 

2.75 

-5.53 

3.15 

HD 227245 


9.81* 

10.25 

9.69 

8.21 

8.07 

7.99 

3.26 

0.91 

2.15 

-4.94 

2.29 

HD 227757 


8.72" 

9.41 

9.27 

8.77 

8.78 

8.76 

3.06 

0.53 

2.15 

-4.02 

2.36 

Cyg OB 8 

HD 191423 


7.43" 

8.19 

8.03 

7.73 

7.72 

7.79 

3.00 

0.46 

2.90 

-5.66 

2.82 

HD 191978 


7.38" 

8.15 

8.03 

7.76 

7.80 

7.80 

2.70 

0.48 

2.75 

-5.46 

2.84 

HD 193117 


8.87" 

9.29 

8.77 

7.24 

7.09 

6.97 

3.18 

0.87 

2.25 

-5.76 

2.37 

CygOB9 

HD 194334 


9.38" 

9.54 

8.82 

6.99 

6.82 

6.74 

2.94 

1.13 

1.19 

-4.88 

1.04 

HD 194649 


9.82* 

9.93 

9.07 

6.86 

6.63 

6.51 

3.10 

1.12 

- 

-4.80 

1.20 

HD 195213 


9.29* 

9.55 

8.82 

6.67 

6.44 

6.26 

3.27 

1.13 

1.00 

-4.88 

1.20 


References, a: iMafz-A pellaniz et al. (2004); b: IWesselius et alJd!982t> : c: lGuarinosldl99lh : <i: lFerniel d 1 9831) : e dKrelowski & Strobell(I2012T) . 


HD 193595: The fit of this object is of good quality. The core 
of the Her line is filled in with some emission that is likely to be 
partially of nebular origin. The stellar parameters agree with a 
main-sequence luminosity class as established in Paper I. 

HD 193682: We use FUSE and IUE spectra to constrain the 
wind parameters (see Table [3]). The FUSE spectrum is correctly 
fitted by our CMFGEN model. The IUE spectrum was observed 
with a low dispersion and the small aperture. Therefore, the qual¬ 
ity of the observation is very low and we must strongly convolve 
the synthetic spectrum to achieve such a quality. The Civ 1548- 
1550 doublet appears too saturated in comparison to the observa¬ 
tion, whilst the fit of the Niv 1718 line seems acceptable. Finally, 
the fit in the optical range is of excellent quality, even though the 
He 1 5876 line is too strong compared to the observations. Unlike 
iRenolust et al.l d2004i) . the red wings of the He n 4686 and Ha 
lines are rather well reproduced. The line profiles show an im¬ 
portant macroturbulence. We find that a combination of v sin i 
of about 150 km s _1 with a macroturbulent velocity of about 
70 km s _1 leads to much better fits than do purely rotationally 
broadened profiles. 

HD 194094: No UV spectrum is available for this star, we 
therefore restrict ourselves to the optical range. The core of the 
Balmer lines is too deep in our synthetic spectra except for Ha. 
Moreover, the He n 4686 is poorly fitted and appears stronger in 
the synthetic model. 

HD 228841: A good fit of the optical wavelength domain is 
achieved, except for the He n 4686 and the He i 5876 lines. The 
best fit is obtained with profiles only broadened by rotation, i.e. 
without any macroturbulence. 

HD 194280: The Aurelie spectrum only allows us to deter¬ 
mine the stellar parameters. The spectral classification is sim- 
ilar to an OC9.7I star according to the tables of iMartins et alJ 
(12005al) . From the CMFGEN best-fit model, we find a C over¬ 
abundance and a N depletion, estimated to about 3.3 and 


0.2 times the solar abunda ncefl respectively, as reported by 
IWalborn & Howarthl ( 2000 ). To obtain a fit of good quality, we 
also need to increase the helium abundance to (He/H) = 0.18 
(by number) i.e. almost twice the solar value ((He/H) 0 = 0.10). 
However, additional diagnostic lines are required to confirm 
such abundances. 

4.2. The Cyg OB3 association 

HD 190864: A good quality of fit is achieved. We need a vol¬ 
ume filling factor of 0.06 and a/? of 1.5 to reach a good accuracy 
for the wind parameters in both UV and optical domains. As for 
HD 193514, Mw„ l -JJ i s consistent with the value determined 
bv lRepolust et alJ (120041) . Moreover, as for HD 193514, the syn¬ 
thetic model fails to correctly fit the N m 4634-41 lines. 

HD 227018: The [4450 - 4900] A wavelength domain is 
properly fitted, which gives us a very good determination of the 
stellar parameters. 

HD 227245: The synthetic model correctly reproduces the 
observations which leads to a good accuracy on the determina¬ 
tion of the stellar parameters. Furthermore, the Hen 4686 and 
Ha lines are well fitted, thereby corresponding to a reliable de¬ 
termination of the mass-loss rate of the star on the basis of these 
diagnostic lines.. 

HD 227757: We need a slow rotation rate to correctly re¬ 
produce the photospheric lines. We also see that the On lines 
are present in the spectrum of the star. The stellar parameters, 
listed in Table [3] are in agreement with those determined by 
IMartins et alJ J2014I) . The luminosity derived for that object and 
its mass-loss rate determined from the Ha line could indicate 
that HD 227757 is a w eak-wind star (see IMartins et al.l l2005bl 
and Martins et ah 2009 for other examples of weak-wind O-type 

1 (C/H) 0 - 2 .4 5 1(U* and (N/H) 0 = 6.02 lO^ 5 (expressed by num¬ 
ber. (Grevesse|etai22007j) 
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stars). However, without any UV spectrum for this star, we can¬ 
not confirm this assumption. 

4.3. The CygOB8 association 

HD 191423: We need to broaden the synthetic model by a 
rotational velocity of 410 km s 'to reach a sufficient quality of 
the fit. Moreover, the model also yields, through the adjustment 
of the carbon lines at 4070 and 4650 A, a large depletion in C 
(about 0.1 times the solar abundance, i.e. 3.5 x 10 s in num¬ 
ber) and an enrichment in N (about 9 times the solar abundance, 
i.e. 5.4 x 10 4 in number). The wind parameters derived on 
the UV P-Cygni profiles indicate a low terminal velocity (about 
600 km s' 1 ). We also see two weak emissions in the external 
wings of Ha. These emissions could be created by a disk of mat¬ 
ter ejected by the high rotational velocity of the star. We stress 
that the P v 1118-1128 and N iv 1718 lines appear slightly weak 
in our model in comparison to the observed spectra. 

HD 191978: The fit is of an excellent quality, except for the 
red wing of the Civ 1548-50 lines and for the Hen 4686 line 
which is too weak compared to the observations. As we already 
stressed, the Ov 1371 and Niv 1718 lines are not good diag¬ 
nostic lines to estimate the clumping factor in late O-type stars. 
Furthermore, given the lack of FUSE spectrum for this star, 
we adopt / = 1.0 (homogeneous model) for HD 191978. We 
need a rotational velocity similar to the macroturbulent velocity 
(v sin i = 40 km s' 1 and v mac = 32 km s' 1 ) to correctly reproduce 
the line profiles. 

HD 193117: The cores of the Balmer lines are too deep in 
the best-fit model. This could come from a nebular emission at 
least for Ha. The other lines are rather well fitted. We also stress 
that although this star appears to be evolved, its nitrogen content 
is solar. 

4.4. The Cyg OB9 association 

HD 194649: In Paper I we determined a brightness ratio of 
4.7 ± 0.4. The fit of the primary disentangled spectrum is good, 
except for the core of the Hfl line. The secondary component be¬ 
ing fainter in the composite spectra, its spectrum appears much 
noisier when corrected for the relative brightness and the disen¬ 
tangling process fails to correctly reconstruct the wings of H/j. 
Nonetheless, the lines are relatively well reproduced by the syn¬ 
thetic model. Although the spectral classifications indicate an 
evolved primary star and a main-sequence secondary star, log g 
determined on the basis of disentangled spectra fail, once again, 
to characterize this luminosity class. 

HD 194334: A good fit is achieved for this object. Only the 
cores of the Balmer lines and N hi 4634-41 are less well repro¬ 
duced. The fits of the Hen 4686 and Ha lines indicate a good 
determination of the mass-loss rate of the star. 

HD 195213: As for several stars in our sample, the cores of 
the Balmer lines are poorly reproduced. However, Ha shows an 
emission component in its centre which probably originates from 
nebular emission. Furthermore, to reach a good fit of the helium 
lines, we have to increase the helium abundance to 0.15 in num¬ 
ber. 


5. Discussion 

5.1. Distances and ages of the OB associations 

The stellar parameters derived in Sect. [4] provide the positions 
of the stars in the Hertzsprung-Russell (HR) diagram. However, 


these locations are dependent on the real distances of the stars. 
As we have already mentioned, the luminosities were com¬ 
puted in the present paper from the distance modulus given by 
[Humphrey s 8 ( 1197 8) and refined through the fit of the SED of each 
object. When no initial value of the distance of the stars is avail¬ 
able in the Humphreys catalogue, we have taken the mean value 
of the correspo nding OB association. Th e observed V and (B-V) 
values are from lKharchenko et al. (j20 09l). whilstHhe bolometric 
corrections were taken from Martins & Plezl (120061) as a func¬ 
tion of the spectral classification of the stars. The error-bars on 
the luminosities of the presumably single stars are mainly de¬ 
termined by the uncertainties on the distances. For the binaries, 
an additional uncertainty is also provided by the brightness ratio 
between both components of the system. 


For this analysis, we use the evolutionary tracks of 
[Brott et all (1201 lab computed for single stars. Given that rota¬ 
tion plays a major role in the evolution of the stars, we select 
the tracks computed with initial rotational velocities of about 
100 km s' 1 (black lines in Fig. Q} and of about 400 km s' 1 
(red lines in Fig. Q3. We stress that an initial rotational veloc¬ 
ity of 300 km s' 1 is generally considered as a ’’standard” value 
for O-type stars. In this analysis, we do not investigate all the 
sample stars together because they are not from the same form¬ 
ing regions and so do not constitute a homogeneous O-type star 
population. 


The distances taken from Humphreys (11978b are close to 
the largest estimates of lUvamker et al. ( 200lb . This yields the 
largest values to the luminosities of the stars. By comparing the 
luminosities computed in the present paper with the observa¬ 
tional values (considered as standard) reported by [Martins et al.1 
(l2005ah . we see that the stars belonging to CygOB9 as well 
as HD 227245 or both components of HD 228989 are fainter 
than the standard expectations. On the contrary, HD 193117 and 
HD 191423 belonging to the CygOB8 association are brighter 
than expected. This makes the estimation of their distance ques¬ 
tionable. It is clear that these standard values are dependent on 
other stellar parameters that are not the same for the stars in our 
sample such as, e.g., logg. To match the standard values for the 
luminosity, we should increase the mean distance of Cyg OB9 to 
a value of 2.0 kpc and decrease that of Cyg OB8 to a mean value 
of 1.9 kpc. 


These assumptions could have a great impact on the loca¬ 
tion of the stars in the HR diagram and could thus affect the 
real determination of their age. Indeed, from these possible dis¬ 
tances for CygOB8 and CygOB9 (about 1.9 and 2.0 kpc, re¬ 
spectively), the stars of CygOB8 could be slightly older (be¬ 
tween 4 and 6 Myrs) than what we find in Fig. |l(c)| whilst those 
of Cyg OB9 could have an estimated age between 2 and 4 Myrs, 
i.e. slightly younger. Furthermore, the location of the secondary 
star of HD, 194649 would be more realistic than what we see in 
Fig. |l(d)| This means that the ages of the stars in Cyg OB8 do 
not seem to confirm the e stimate s of ^_3 Myrs computed from 
the evolutionary tracks bv llJvaniker et al.l (12001 ) for this associ¬ 
ation, but this assumption must still be verified on the basis of a 
larger sample of stars belonging to this association. 

From the distances quoted in Table [2] we estimate the evo¬ 
lutionary ages of O-type stars in Cyg OBI and in CygOB3 be¬ 
tween 2 and 7 Myrs and between 3 and 5 Myrs, respectively, ac¬ 
cording to the parameters reported in Tabled We must be careful 
about these HR diagrams, because the real ages of the stars can 
be slightly different depending on their initial rotational velocity. 
Therefore, some stars such as HD 193682 or HD 228841 (with 
vsin; = 150 and 317 km s' 1 , respectively) could be older than 
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Table 3. Stellar and wind parameters derived for the O-type stars in four Cygnus OB associations. 


Source 

SpT 

Teff 

lo s& 

logg 

loggtme" 

R 

log {M 

,heo ) 6 

io «%>r 

U 

Voo 

P 

v sin i 

^mac 

M evo i^ 

M S p ec c 


N/ 



[kK] 




[Re] 




[km s~'] 


[km s~*] 

[km s~'] 

[M g ] 

[M G ] 

[Ms} 

[Xl0“‘ 










CygOBl 








P 

cr 


HD 193443P 

0 9111 

31.2 

5.35 ±0.13 

3.65 

3.66 

16-2^43 

~ 

- 6.8 

- 

- 

- 

- 

105 

38 

30.04 3 

43 8 +1U3 U 

-to.°_30 4 

1-0 ±04 

0.6 ± 0 . 

HD 193443S 

0 9.5 V 

29.7 

4.76 ± 0.23 

3.75 

3.76 

9-1-13 


- 

- 

- 

- 

- 

100 

27 

18-7^53 

17 4 + 57.2 

1 ' -13.2 

0.4 ± OH. 

0.6 ± 0 . 

HD 193514 

O 7-7.5111(f) 

34.5 

5.65 ± 0.09 

3.50 

3.51 

18.7-1:1 


-6.5 

-5.6 

0.03 

2190 

1.5 

90 

27 

45.7*3 

41-31 173 

=r 

2.0 ± 1 . 

HD 193595 

0 7V 

35.5 

5.11 ±0.05 

3.75 

3.75 

9 5 +L2 


- 

-8.7 

1.0 

- 

- 

50 

28 

26.5+23 

18- 7 -6°8 9 

c 1 

3.0 ± 1. 

HD 193682 

0 5111(f) 

39.4 

5.50 ± 0.09 

3.75 

3.77 

12 1 +20 
az " a -1.7 


-5.8 

-5.7 

0.01 

2650 

1.0 

150 

73 

4 1 -4+4j 

31 4+22.5 
:7A -13.1 

& 

p 

6.0 ±3. 

HD 194094 

0 8 III 

29.9 

5.45 ± 0.36 

3.25 

3.26 

19-8!‘ 2 6 3 


-6.7 

-7.0 

1.0 

- 

- 

65 

32 

33 o +23 - 8 

25 9 +59 - 7 
-18.0 

3 

or 

1.0 ±0. 

HD 194280 

0 9.71 

26.9 

5.61 ±0.36 

3.10 

3.11 

9Q 4+18.6 
z '' 7 ' -11.4 

< 

-6.5 

- 

- 

- 

- 

92 

41 

oq 9+25.4 
-12.9 

41 1+96.8 
' -28.8 

£ 

0.1 ±0. 

HD 228841 

07 

34.5 

5.42 ± 0.37 

3.50 

3.62 

14.4+^;“ 

~ 

-6.4 

-7.1 

1.0 

- 

- 

317 

0 

99 '7+22.7 
^ -10.1 

31 4 +72 - 9 
JA -21.9 

p- 

5.4 ±2. 

HD 228989P 

0 8.5V 

29.4 

4.43 ±0.14 

3.50 

3.56 

6 3 +26 


- 

- 

- 

- 

- 

145 

3 

15.91 3 ; 2 

5 3 +13 - 3 

7.0 ± 0 | 

1.0 ±0. 

HD 228989S 

0 9.7 V 

27.8 

4.28 ±0.15 

3.75 

3.78 

6 0 +26 

U -1.7 


- 

- 

- 

- 

- 

120 

3 

13-0+ 3 ;4 

7 7+20.6 
'■'-5.6 

6.1 ± 0 l 

C/D 

3.0 ± 1. 

HD 229234 

0 9III + ... 

32.0 

5.09 ± 0.20 

3.60 

3.61 

11 4 +3 - 9 

11 •‘+_2 9 

~ 

- 6.8 

-6.5 

1.0 

- 

- 

85 

38 

24-8+43 

1 Q 4+24.6 
A ^-10.8 

O 

b 

r\ 

0.6 ± 0 . 










Cyg OB 3 










HD 190864 

0 6.5111(f) 

38.0 

5.35 ±0.14 

3.75 

3.76 

10.9^ 


-6.3 

-6.4 

0.06 

2250 

1.0 

80 

43 

34 1 + * A 

94 Q+23.2 

^•^-12.0 

-d- 

O) 

5.0 ±3. 

HD 227018 

0 6.5 

37.3 

5.50 ±0.18 

4.00 

4.00 

13.5*° 


-6.9 

- 

1.0 

- 

- 

70 

35 

77 9+13.4 

J ' * -4.8 

66 - 7 !ii 

C/D 

1.0 ± 0 . 

HD 227245 

0 7111 

36.0 

5.18 ±0.14 

3.75 

3.75 

10 -O^jj 

~ 

-6.4 

-6.7 

1.0 

- 

- 

65 

19 

21J -23 

20 . 8 :|” 

■ 8 - 

3.0 ±0. 

HD 227757 

0 9V 

32.0 

4.73 ±0.13 

3.75 

3.75 

n /C+ 1.8 

7 ' 6 -1.4 


- 

- 8.2 

1.0 

- 

- 

45 

3 

18.6^0 

11 . 8!‘ 0 6 9 

P, 

P 

0.6 ± 0 . 










Cyg OB 8 








o' 


HD 191423 

ON 9IIIn 

30.6 

5.42 ± 0.28 

3.50 

3.67 

18.3-“;' 

~ 

- 6.8 

-6.4 

0.01 

600 

0.9 

410 

0 

34J-- 

54 6 +9i l 

-34.6 

5* 

5.4 ±2. 

HD 191978 

0 8 III 

33.2 

5.35 ± 0.23 

3.75 

3.75 

14 3 +5 - 5 

A ^ -4.0 


-6.7 

-8.7 

1.0 

1600 

- 

40 

32 

79 1 +9.4 
a -7.4 

49 9+59.8 
^ z " z '-24.7 

O? 

C 

3.0 ± 1. 

HD 193117 

0 9III 

30.1 

5.36 ±0.11 

3.30 

3.31 

17.6*;* 

~ 

- 6.8 

-6.3 

1.0 

- 

- 

80 

28 

30.4! 3 : ; 

97 9+19.9 
-10.7 

n 

0.6 ± 0 . 










Cyg OB 9 








OQ 

3 


HD 194334 

O 7-7.5111(f) 

33.9 

5.19 ±0.07 

3.50 

3.51 

11 4 +L/ 

~ 

-6.5 

- 6.2 

1.0 

- 

- 

62 

32 

26.61+;: 

15 3 +lu - 2 

AJ -6.1 

c»- 

3.0 ± 1. 

HD 194649P 

0 6111(f) 

38.0 

5.14 ± 0.15 

3.75 

3.75 

8 . 6 +H 


- 6.1 

- 

- 

- 

- 

60 

67 

28.1*3 

15 3 +35 1 
AJ -10.6 

4.9 ± o§ 

3.0 ± 1. 

HD 194649S 

0 8 V 

33.0 

4.36 ± 0.27 

3.75 

3.76 

4.6+jl 


- 

- 

- 

- 

- 

60 

27 

9 6 +9A 

y ‘ KJ - 9.0 

4 5 +15 - 9 

^ -3.5 

1.9 ±0l 

O 

1.0 ± 0 . 

HD 195213 

0 7111(f) 

35.3 

5.19 ±0.22 

3.50 

3.51 

i°-5il:l 

~ 

-6.4 

- 6.0 

1.0 

- 

- 

85 

28 

27 6 +5 - 7 
^ /,u -6.1 

1 9 9+17.8 

6 

2. 

p 1 

6.0 ±3. 


o 


Notes. ,a> g true : gravity corrected from the contribution of the centrifugal force i.e. g true = g + (v sin i) 2 /R (b) log(M t hco: theoretical mass-loss rate estimated bv lMuiires et alj j2012h for a given spectral “ 

classification (c) \og(Mn a \ mass-loss rate (where the correction of the clumping yff is not applied) computed from the H„ line. (rf ' M eV oi: mass estimated from the HR diagram (e) M spec : mass 
computed from M = gR 2 /G <2> M dyn : minimum mass of a binary component (M 3 sin i) 
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Fig. 1. Positions of the Cygnus O-type stars in the HR diagram. Evolutionary tracks and isochrones are from iBrott et ali ( 201 lab 
and were computed for v rotjni[ = 100 km s^ 1 (black lines) and for v ro t inil = 400 km s' 1 (red lines). Isochrones, represented by dashed 
lines, correspond to ages ranging from 2 to 8 Myrs with a step of 2 Myrs. 


inferre d from their position in the HR diagram (see iBrott et al] 
1201 laL for further details). 

Determining the real positions of the O stars from the HR 
diagram is difficult and requires a good knowledge of their real 
distances. Therefore, our study stresses the importance of astro- 
metrical missions such as Gaia to better understand the evolu¬ 
tionary properties of massive stars. 

5.2. The N content 

Nitrogen is an important element to analyse and to understand 
the rotational mixing in the evol utionary processes of a mas¬ 
sive star. Theoretical studies (e.g. lMevnet & MaedeifeOOd) have 
shown that rotationally induced mixing can affect the chemical 
composition of the surface layer of massive stars even though 
they are still on the main-sequence band. The theory thus pre¬ 
dicts that the larger the rotational velocity of a star, the higher 
its nitrogen abundance. Several observational studies on B stars 
in the LMC (notably the ESO VLT-FLAMES survev lEvans et al.1 
[20051 . 2006 1 have shown that at least 60% of the B stars dis¬ 
play the pattern expected from evolutionary models which in¬ 
clude the rotational mixing. However, simulations of LMC early 


B-type stars made by iBrott et al. (201 lb l failed to reproduce 
the remaining 40%, which are split into two different groups 
of stars, one containing slowly rotating, nitrogen-enriched ob¬ 
jects and another one containing rapidly rotating un-enriched 
objects. The analysis of iRivero Gonzalez et al.l (1201 2l) also in¬ 
dicated a large number of O-type stars in the LMC with large 
enrichment in nitrogen and with a low rotation rate, support¬ 
ing the same conclusions as for the B stars in the LMC. These 
authors stressed, however, that the problem of the group con¬ 
taining the highly-enriched slow-rotating stars was more severe 
for the O-type star population than for the B-type star one. A 
possible explanation concerning the presence of the group with 
high-enriched slow-rotating stars was advanced bv IMorel et al.l 
(i2006t) who have analysed ten slowly rotating galactic early [3- 
Cephei B-type stars and found that out of four heavily enriched 
stars, three had a magnetic field. 

In the present paper, we determine the N content of nine¬ 
teen O stars belonging to four Cygnus OB associations. This 
result is combined and put in perspective with several other in¬ 
vestigations of O-t ype stars in the Galaxy dMartins et al.ll2012l 
iBouret et al.ll2012l and lMartins et al.ll2014l) . Although this con¬ 
stitutes a set of heterogeneous studies, all the analyses use 
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the same tool to determine the stellar parameters. We show in 
Fig. El the nitrogen content as a function of the projected ro¬ 
tational velocities of 90 O-type stars, the so-called Hunter dia¬ 
gram dHunter et alJ 12009 1. We stress that we remove from our 
sample all the binary components as well as a few st ars listed in 
Marti ns et al.l (120141) that are binaries, e.g. HD 93250 (ISana et al.1 
2011 h or HD 193443 (Paper I). We tentatively define possibl e 
groups similar to those introduced by [Hunter et al.l ( 82008112009 1. 
To do that, the group constituted of stars with intermediate rota¬ 
tion rates (between 50 and 110 km s 1 j and a high enrichment in 
nitrogen is assumed to be empty as noted by [Brott et al.l d2011 bl) . 
In this diagram, there is a clear outlier (HD 194280) that seems 
to belong to none of the different groups. We stress though that 
defining their exact location requires a_more sophisticated the¬ 
oretical work similar to that of [Brott et al.l d201 lbl) for massive 
stars in the LMC. 



Fig. 2. Nitrogen surface abundance (in units of 12 + log(N/H)) as 
a function of projected rotational velocity, the so-called Hunter 
diagram. O stars in our sam ple are shown in black, the O- 
sta rs of NGC 2244 (IMartins et alJ l2012h in red, stars studied 
byjBouretetal:J 2012h in blue, and the MiMes O-type stars 
dMartins et all2014l) in green. Grey lines define the possible lo¬ 
cations of the different groups. 


In Fig. El we show the N content of our initial targets ver¬ 
sus their luminosity. We clearly see a trend of higher N_content 
forjnore luminous stars. As already stressed by [Martins et all 
(2012), the fast rotators do not appear as outliers. Two clear out¬ 
liers are shown in this figure: HD 194280, as we have already 
mentioned, with a very low nitrogen abundance and, to a lesser 
extent, the secondary component of HD 228989. 

5.3. Stellar masses 

The differences between the evolutionary masses, i.e. the masses 
determined from the HR diagram and the evolutionary tracks, 
and the spectroscopic masses, i.e. those determined from the 
atmospherical parameters (7 d r, log#, and L), are recurrent is¬ 
sues in astrophysics. This so-called ’’mass discrepancy” refers 
t o a sys t emat ic overestimate of the former relative to the latter. 
iRepolust et al.l ( 2004 ) already noted that for stars with masses 
smaller than 50 M 0 a parallel relation to the 1:1 relation in 
the M evo i - M s p e c diagram could be followed. More recently, 
iMartins et al . ('2012) . from a small sample of O-type stars in 



Fig. 3. Nitrogen surface abundance (in units of 12 + log(N/H)) 
as a function of stellar luminosity. Dots represent the single 
O stars in our sample, whilst triangles are the binary compo¬ 
nents. Evolutionary tracks are (from left to right) computed for 
stars with 10, 15, 20, and 40 M 0 and initial velocities of about 
300 km s -1 . 


NGC 2244 and in Mon OB2, reported on a clear trend of mass 
discrepancy for stars with M < 25 M 0 , consistent with the re¬ 
sults of IRepolust et al.l d2004l) . 

In Paper I we also determined the dynamical mass of each 
component of the binary systems. We present them in Table E] 
as a reminder. We can thus compare them to their evolutionary 
and spectroscopic masses. The evolutionary masses have been 
estimated for each star (presumably single or individual compo¬ 
nent) by a bi-interpolation of the evolutionary tracks on T e ff and 
on log (L/Lq), whilst the spectroscopic masses have been com¬ 
puted from M = gR 2 /G, where G is the universal gravitational 
constant. We stress that these masses were computed from the ef¬ 
fective gravities corrected for the effects of the centrifugal forces 
caused by rotatio n (gn-ii p in Table E} by following the approach 
of [Herrero et al.l d!992l) . We clearly see a real disagreement for 
the primary star of HD 228989. First, the spectroscopic mass is 
smaller than the dynamical mass that is supposed to be the min¬ 
imum mass. Then, the mass ratio computed from the spectro¬ 
scopic values indicates that the secondary component is more 
massive than the primary one which is the opposite of the results 
provided in Paper I. This assumes that the distance and/or the 
gravities of both components are poorly estimated. The reason 
for a poor determination of the log g is that the disentangling pro¬ 
cess does not perfectly reproduce the wings of the Balmer lines 
because of the broadness of these lines and because the two in¬ 
dividual line profiles are not completely resolved throughout the 
orbit. 

When we only focus on the single stars, we find that the ob¬ 
jects in our sample (Fig. [4} with a mass smaller than 30 M 0 
present a clear mass discrepancy as already detected for O- 
type s tars belonging to NGC 2244 and Mon OB2 (IMartins et al.l 
|2012j)- For the stars in our sample with a higher mass, the un¬ 
certainties do not allow us to conclude whether or not the mass 
estimates are significantly different. Two ingredients play an im¬ 
portant role in the determination of these uncertainties: the dis¬ 
tance and log#. As we have already mentioned, the distances of 
the stars are still poorly known. It affects the stellar luminosities 
and thus the stellar radii since the T e s is thought to be well con- 
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Fig. 4. Comparison of the evolutionary masses and the spectro¬ 
scopic masses of the presumably single stars. 


strained. Concerning log#, we made several tests to correctly fit 
the wings of the Balmer lines on echelle spectra and also on the 
Aurelie spectra (see Paper I, for a complete list of the data). Our 
estimation of the stellar gravity thus seems quite robust, there¬ 
fore minimizing the contribution of the log g to the uncertainty 
budget of M spec . 

5.4. Mass-loss rates 

We also estimate the mass-loss rates (Table[3j for most objects in 
our sample, but the lack of UV spectra does not allow us to deter¬ 
mine all the wind parameters of these stars. Often, the terminal 
velocity, /, and the /j parameters cannot be estimated from only 
the Her line. Moreover, when no FUSE data were available, we 
only assumed that the stellar winds are homogeneous (/ = 1.0). 

The mass-loss rates obtained from atmosphere models and 
included in Table [3] can be compared to the theoretical val¬ 
ues quoted by iMuiires et al.l ( 201A These values are given as 
a function of spectral type, but some of them are missing no¬ 
tably for late main-sequence stars. For most of our objects for 
which a homogeneous wind was assumed, the agreement be¬ 
tween estimated and theoretical values is relatively good. Only 
few stars (HD 228841, HD 193514, HD 227018, HD 193117, 
and HD 191978) have mass-loss rates which differ by more than 
about 0.4 dex. Although the quantities M/ yff remain quite close 
between the observations and the theory, we can see that the ob¬ 
served mass-loss rates in general appear slightly larger than the 
theoretical values. 

6. Conclusion 

We analysed by means of the CMFGEN atmosphere code, the 
fundamental properties of fifteen presumably single stars and 
of the individual components of four binary systems belong¬ 
ing to different Cygnus OB associations: CygOBl, CygOB3, 
CygOB8, and CygOB9. In addition to the optical spectra, we 
retrieved, for several single stars, the UV spectra (FUSE and/or 
IUE) to better constrain the wind parameters of these objects. 

The analysis of the individual parameters showed that the 
luminosities of several objects do not agree with the standard 
luminosities of stars with similar spectral classification. This is 
most likely because the distances of these stars are poorly con¬ 


strained. Therefore, it is more likely that the positions of cer¬ 
tain stars should be shifted upwards in the HR diagram. A direct 
consequence is that the stars in Cyg OB8 and Cyg OB9 could be 
younger than we have determined. However, if the distance used 
in our analysis is correct, this wouldjnean that Cyg OB8 is older 
than 3 Myrs as quoted bv lUvamker et alj d2001l) . 

We compared the N content of the stars in our sample 
to other galactic O-type stars. The resulting Hunter diagram 
shows the same structure (the sam e five groups) as f or the 
O and B-type stars in the LMC dHunter et al.l 120081 120091 : 
iRivero Gonzalez et alJl2012h . 

By comparing the evolutionary and the spectroscopic 
masses, we also pointed out a discrepancy between the masses 
for the stars in our sample. This overestimate of the evolutionary 
masses relative to the spectroscopic masses is particularly obvi¬ 
ous for stars with mas ses smaller than 30 M Q , as was already 
noted by iMartins et al.1 d2012l) . Above this value, the masses 
agree with each other within the uncertainties. Although this 
sample of stars takes into account objects_more evolved than 
those in the investigation of lMartins et alj d2012l ). we obtain sim¬ 
ilar general conclusions. 
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Fig. 5. Synthetic SEDs (in red) compared to UBVJHK photometry for the stars in our sample (in black). The distance, E(B-V) and 
Ry were determined from the best agreement between models and observations. 
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Fig. 6. (continued). 
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Fig. 7. Best-fit model for HD 193443 primary and secondary (red line) compared to disentangled spectra (black line). 















Normalized flux 


L. Mahy et al.: The fundamental parameters of O-type star population in four Cygnus OB associations. Online Material p 4 


~\ -1-r 


t -1-r 


t -1-r 


t -1-r 


2 - 


1 


J_I_l_ 


J_I_L_ 


J_I_L 


J_I_l_ 



1250 


1300 


1350 


1400 


1450 



t -r 


t -1-1-r 


t -1-1-r 


2 - 


1 


J_L 


J 

1650 


J_I_I_L 


J_I_I_L 


1700 

A [A] 


1750 


Fig. 8. Best-fit model for HD 193514 (red line) compared to FUSE and IUE spectra (black line). 
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Fig. 9. Best-fit model for HD 193514 (red line) compared to Elodie spectrum (black line). 
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Fig. 10. Best-fit model for HD 193595 (red line) compared to Espresso spectrum (black line). 
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Fig. 11. Best-fit model for HD 193682 (red line) compared to FUSE and IUE spectra (black line) 
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Fig. 12. Best-fit model for HD 193682 (red line) compared to Espresso spectrum (black line) 
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Fig. 13. Best-fit model for HD 194094 (red line) compared to Espresso spectrum (black line). 
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Fig. 14. Best-fit model for HD 194280 (red line) compared to Aurelie spectrum (black line). 
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Fig. 15. Best-fit model for HD 228841 (red line) compared to Espresso spectrum (black line). 
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Fig. 16. Best-fit model for HD 228989 primary and secondary (red line) compared to disentangled spectra (black line). 
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Fig. 17. Best-fit model for HD 229234 (red line) compared to Espresso spectrum (black line). 
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Fig. 18. Best-fit model for HD 190864 (red line) compared to FUSE and IUE spectra (black line). 
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Fig. 19. Best-fit model for HD 190864 (red line) compared to Elodie spectrum (black line). 
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Fig. 20. Best-fit model for HD 227018 (red line) compared to Aurelie spectrum (black line). 
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Fig. 21. Best-fit model for HD 227245 (red line) compared to Espresso spectrum (black line). 
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Fig. 22. Best-fit model for HD 227757 (red line) compared to Espresso spectrum (black line). 
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Fig. 23. Best-fit model for HD 191423 (red line) compared to FUSE and IUE spectra (black line). 
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Fig. 24. Best-fit model for HD 191423 (red line) compared to Espresso spectrum (black line). 
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Fig. 25. Best-fit model for HD 191978 (red line) compared to IUE spectrum (black line). 
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Fig. 26. Best-fit model for HD 191978 (red line) compared to Elodie spectrum (black line). 






















Normalized flux 


L. Mahy et al.: The fundamental parameters of O-type star population in four Cygnus OB associations. Online Materialp 23 


1.1 
1 

0.9 
0.8 
0.7 
0.6 

4000 4100 4200 4300 4400 

1.1 
1 

0.9 
0.8 
0.7 
0.6 

4500 4600 4700 4800 4900 

1.05 
1 

0.95 
0.9 
0.85 

5350 5400 5450 ~ 6520 6540 6560 6580 6600 

1.1 

1 

0.9 
0.8 
0.7 

5600 5700 5800 5900 

xrAi 

Fig. 27. Best-fit model for HD 193117 (red line) compared to Espresso spectrum (black line). 
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Fig. 28. Best-fit model for HD 194334 (red line) compared to Espresso spectrum (black line). 
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Fig. 29. Best-fit model for HD 194649 primary and secondary (red line) compared to disentangled spectra (black line). 
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Fig. 30. Best-fit model for HD 195213 (red line) compared to Espresso spectrum (black line). 
































